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In this note we present several recent results for the production of multijet final states in pp 
collisions at a center of mass energy of 1.96 TeV. These measurements, which include the cross- 
section for dijet and three-jet production and the ratio of three-jet to two-jet cross sections, were 
taken with the D0 experiment at the Fermilab Tevatron collider. 

I. INTRODUCTION 

In hadron-hadron collisions, production rates of coUiniated sprays of hadrons, called jets, are sensitive to both 
the dynamics of the fundamental interaction and to the partonic structure of the initial-state hadrons. The 
latter is usually parameterized in parton distribution functions (PDFs) of the hadrons. Thus measurements 
of jet production properties can be used to test the predictions of quantum chromodynamics (QCD) as well 
as to constrain the PDFs. Because jet production has the largest cross-section of any high px process at a 
hadron collider, the study of jet production provides the highest energy reach for the experiment and a unique 
sensitivity to new physics. 

Several jet measurements are now well established at the Tevatron, including the measurement of the inclusive 
jet cross section and of the angular and invariant mass dependencies of dijet events These measurements 

provide important tests of perturbative QCD (pQCD) and are input to PDF calculations. Measurements of 
multi-jet production take advantage of the fact that these processes have the same PDF sensitivity as dijet 
production, but are sensitive to processes to third order in the strong coupling constant as- Fundamental to 
the understanding of order processes is the three-jet cross-section, which we present as a function of the 
invariant mass of the three-jet system. Studies dedicated to the dynamics of the interaction are preferably 
based on observables which are insensitive to the PDFs. Such observables can be constructed as ratios of cross 
sections for which the PDF sensitivity cancels. In this note we report a measurement of -R3/2, the ratio of the 
inclusive three-jet to the inclusive 2-jet cross-sections. 

In this note we present several recent results for the production of multijet final states in pp collisions at a 
center of mass energy of 1.96 TeV. These measurements, which include the cross-section for dijet and three-jet 
production and the ratio of three-jet to two-jet cross sections, where taken with the D0 experiment [l| at the 
Fermilab Tevatron collider. The data analyzed were recorded using a single jet triggers at a variety of thresholds, 
corresponding to an integrated luminosity of 0.7 fb— 1. The event selection, jet reconstruction, jet energy and 
momentum correction in this measurement follow closely those used in our recent measurements of inclusive jet 
and dijet distributions @, 0]- The primary tool for jet detection is the finelysegmented uranium- liquid argon 
' calorimeter that has almost complete solid angle coverage 1.7° <0 < 178.3° mi. Jets are defined by the Run II 
midpoint cone jet algorithm ^] with a cone radius (for most jet studies) of i?cone = a/ (Ay)^ + (A^)^ = 0.7 in 
rapidity y and azimuthal angle (f>. Rapidity is related to the polar scattering angle 9 with respect to the beam 
axis hy y = 0.5 In [(1 -I- /3cos0)/(l — Pcosd)] with /3 — \p\/E. The jets in an event are ordered in descending 
transverse momentum px with respect to the beam axis. 

For all jet studies, a jet energy scale (JES) correction is applied to jet four momentum as the measured 
four momentum of a jet is not the same as that of a jet entering the calorimeter due to the response of the 
calorimeter; energy showering in and out of the cone; and additional energy from detector noise, event pile-up 
and multiple pp interactions. The JES correction is determined using the pt imbalance in 7-f jet and dijet 
events. The additional energy from pile-up and multiple interaction is determined from a minimum bias sample. 
The JES corrections are of the order of 50% for a jet energy of 50 GeV and 20% for a jet energy of 400 GeV. 
However, for the 0.7 pb— 1 data sample used for the measurements presented in this talk, collected with the 
D0 detector during 2004-2005 in Run II of the Fermilab Tevatron Collider, the uncertainty on the JES is very 
small, less than 2% over jet p^'s from 60 to 300 GeV. 

For the sake of brevity, the reader is directed to reference fsl for a discussion of the inclusive jet cross section 
measurement, and reference [3| for the details on the extraction of the strong coupling constant from the inclusive 
jet cross section. 
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FIG. 1: The fractional uncertainty from the jet energy scale (JES) for jets of cone size Rconc ~ 0.7 and central pseudo- 
rapidity. 

II. DIJET CROSS-SECTIONS 

Of the variables that fully specify the kinematics of dijet production in hadron collisions, the dijet invariant 
mass Mjj and rapidity difference y* (or equivalently Xdijet) probe the features of LO 2^2 processes, while 
measurements of the azimuthal opening angle Acj) and ratio of transverse momenta pt2/pti can be used to 
determine the effect of higher-order processes. The dijet invariant mass cross section is particularly sensitive 
to the gluon contribution to PDFs at high x. In addition, the dijet mass spectrum can be used to search for 
evidence of quark compositeness or new strongly produced resonances. 

D0 has measured the dijet cross section as a function of the dijet invariant mass Mjj Q . This cross section 
is measured in six regions of ymax, the rapidity of the leading jet in the event. 

Figure [5] shows the resulting measurements, while figure [5] shows the measured cross section divided by the 
NLO prediction from the NLQJET-i-i- [1] program and using the MSTW2008NLO PDFs 0. This measurement 
in an important extension of pQCD to the forward region, and includes measurements of dijet invariant masses 
up to 1.2 TeV. 
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FIG. 2: Dijet mass cross sections, shown for six regions of Ij/maxI- (Note that these regions scaled for plotting, wtih the 
exception of \y,nax\ < 0.4.) Full lines correspond to the NLO calculations with NLOJET-f-f and MSTW2008 PDFs. 
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FIG. 3: Dijet mass cross section divided by NLO prediction, with total systematic uncertainties indicated by the shaded 
bands. The dataset corresponds to an integrated luminosity of 0.7 fb~^. 



III. THREE JET CROSS-SECTION 

D0 has measured [l^l the differential inclusive three-jet cross section as a function of the invariant three-jet 
mass for three hard, well separated jets in three regions of jet rapidities {\y\ < 0.8, \y\ < 1.6, \y\ < 2.4 for 
all three jets) and in three regions of the third jet transverse momentum {pxs > 40 GeV, pxs > 70 GeV and 
Pt3 > 100 GeV) with leading jet pr > 150 GeV. The three-jet mass spectra are shown in Fig. ??. This 
represents the first fully corrected 3-jet mass distribution to be published, and the first comparison of such a 
cross-section to NLO pQCD calcuations. 
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FIG. 4: a) Three-jet mass cross section in regions of jet rapidities. The \y\ < 2.4 region is scaled by a factor of 4 
for readability. Systematic uncertainty is shown by shaded band. Full lines correspond to the NLO calculations with 
NLOJET+-I- and MSTW2008 PDFs. No events are found in the highest Msjet bin in \y\ < 0.8 region, b) Three-jet 
mass cross section in regions of the third jet transverse momenta. The pra > 40 GeV region is scaled by a factor of 2 
for readability. Systematic uncertainty in all three-jet mass bins is shown by shaded band. Full lines correspond to the 
NLO calculations with NL0JET4-+ and MSTW2008 PDFs. In both figures, the dataset corresponds to an integrated 
luminosity of 0.7 fb~^. 



The data are compared to the perturbative QCD NLO predictions calculated in NLOJET- 
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MSTW2008NLO@ PDFs. The NLO prediction is corrected for hadronization and underlying event effects, with 
correction varied in the range from -3% to 6%, obtained from PYTHIA[ii| simulations using tune QW[ll]- The 
comparison of data to theory is shown in Fig. [5] and [6] 
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FIG. 5: Data to theory ratio of the three-jet mass cross section in three regions of jet rapidities. The total systematic 
uncertainty is shown by a shaded band. The PDF uncertainty comes from the 20 MSTW2008NLO eigenvectors. The 
scale uncertainty is determined by varying the scale up and down by a factor of 2. 
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FIG. 6: Data to theory ratio of the three-jet cross section in three regions of the third jet transverse momentum. The 
total systematic uncertainty is shown by a shaded band. The PDF uncertainty comes from the 20 MSTW2008NLO 
eigenvectors. The scale uncertainty is determined by varying the scale up and down by a factor of 2. 



IV. RATIO OF THREE-JET TO DUET CROSS-SECTION, Tig/a 

The ratio of inclusive jet cross sections, R3/2{Pt'^^) — {dcr>3jct / dp^'^^) / {da>2-jct / dp'i^'^^) , is less sensitive to 
experimental and theoretical uncertainties than the individual cross sections, due to cancelations of correlated 
uncertainties. Here R3/2 is measured as a function of the leading jet pr in an event, p'ji^^, in the interval 
(p?;'" + 30 GeV) < pl^^"" < 500 GeV, for p^'" requirements of 50, 70, and 90 GeV. 

The results are displayed in Fig. [71 where the inner error bars represent the statistical uncertainties while 
the total error bars represent the quadratic sums of statistical and systematic uncertainties. The data are 
compared to the predictions from different Monte Carlo event generators. The SHERPA vl.1.3 predictions [l3| . 
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which include the tree-level matrix elements for 2-, 3-, and 4-jet production, are shown as solid lines in Figs.[71 
These were obtained using default settings and MSTW2008LO PDFs and by matching the leading order matrix 
elements for 2-, 3-, and 4-jet production with a parton shower. 

In addition, the data are compared to predictions from the PYTHIA event generator (version 6.422). The 
matrix elements implemented in PYTHIA are only those for 2-jet production. All additional jet emissions are 
produced by a parton shower. There are two different implementations, a virtuality-ordered parton shower and 
apT-ordered one. Both are highly tunable and more than 50 tunes are provided in PYTHIA v6.422. All tunes 
studies here use the CTEQ5L PDFs 

In Fig. [7] the data are compared to PYTHIA tunes which use the pT-ordered parton shower [H, [l^. These 
are tune "Professor pTO" [13] and two extreme tunes from the "Perugia" series of tunes [l^ : the tunes "Perugia 
hard" and "Perugia soft". All of these tunes give very different results for but all predict significantly 

higher ratios i?3/2 than what is seen in the data, even for the softest tune from the Perugia series. 
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FIG. 7: The ratio i?3/2 of trijet and dijet cross sections, measured as a function of the leading jet pt (j>t^^) for different 
p™'" requirements for the other jets. The predictions of SHERPA and pythia (for Tunes A, BW, and DW (top); and for 
three tunes using the pT-ordered parton shower (bottom)) are compared to the data. 

The data are compared to the perturbative QCD NLO predictions calculated in NLOJET-I--H 8] program with 
MSTW2008NLO 3 PDFs. The NLO prediction is corrected for hadronization and underlying event effects, with 
correction varied in the range from -3% to 6%, obtained from PYTHIA JTJ simulations using tune QWpJi]. The 
comparison of data to theory is shown in Fig. |S| 



V. SUMMARY 



The D0 experiment at the Fermilab Tevatron collider has been enormously successful, and continues to pro- 
duce important physics measurements, including tests of the Standard Model and searches for new phenomena. 
Only a small portion of the QCD-related analyzes we have pursued could be shown in this talk. The reader is 
urged to check http : / / www- dO . f nal . go v/re suit s/index . html for the latest results. 
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FIG. 8: The ratio R3/2 of trijet and dijet cross sections, measured as a function of the leading jet pt (pt'"') for different 
p™'" requirements for the other jets. The predictions of next-to- leading order pQCD, calculated with NLOJETH — h, are 
compared to the data. 
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